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Assuming that the lightest supersymmetric particle (LSP) is the lightest neutralino x, we 



present a detailed exploration of neutralino-stau (x — f) coannihilation channels, including 
analytical expressions and numerical results. We also include x coannihilations with the e 
and fx. We evaluate the implications of coannihilations for the cosmological relic density of 
the LSP, which is assumed to be stable, in the constrained minimal supersymmetric extension 
of the Standard Model (CMSSM), in which the soft supersymmetry-breaking parameters are 
universal at the supergravity GUT scale. We evaluate the changes due to coannihilations 
in the region of the MSSM parameter space that is consistent with the cosmological upper 
limit on the relic LSP density. In particular, we find that the upper limit on m~ is increased 
from about 200 GeV to about 600 GeV in the CMSSM, and estimate a qualitatively similar 
increase for gauginos in the general MSSM. 



1 Introduction 



One of the most appealing candidates for the cold dark matter in the Universe is the lightest 
super symmetric particle (LSP). This is stable if the quantum number R = (—\Y B+L+ j s 
conserved [0, as in the minimal supersymmetric extension of the Standard Model (MSSM) 
0, and hence a candidate relic from the Big Bang. Stringent upper limits on the relative 
abundances of anomalous heavy isotopes suggest that the relic LSP should be electrically 
neutral with no strong interactions, so as to ensure that it does not bind to nuclei. Weakly- 
interacting candidates for the LSP, within the MSSM, include the sneutrinos U{ and the light- 
est neutralino x- LEP limits on Z° — > invisible neutral particles suggest that > M z o/2, 
in which case direct searches for dark matter particles along with cosmological constraints, 
remove any sneutrino from consideration^ as the dark matter in the MSSM [|J]. Thus super- 
symmetric dark matter is commonly thought to consist of x neutralino particles [|. 

It is a remarkable feature of x dark matter that its cosmological relic density naturally 
falls in the range allowed by cosmology and preferred by astrophysics in generic domains of 
MSSM parameter space ||. This is in agreement with general arguments that the mass of a 
cold dark matter particle whose relic density is fixed at freeze-out from thermal equilibrium 



should not be more than ~ y/To x Mp ~ 1 TeV, where To ~ 2.73 K is the present cosmic 
microwave background temperature and Mp ~ 1.2 x 10 19 GeV. This general argument 
suggests that such a dark matter particle should be detectable in experiments at the LHC [|J, 
which is also the conclusion reached by studies of the physics reach of the LHC in the 
parameter space of the MSSM ||. 

It has been appreciated for some time that the relic density at freeze-out may be sensitive 
to coannihilation processes involving the LSP and heavier sparticles, % [[H| [11]. The relative 
importance of such coannihilation effects is controlled essentially by the ratio of coannihi- 
lation and annihilation cross sections: a~^/ <xW a~~ and the ratio of number densities, 
which is determined by a Boltzmann factor: exp(m-^ — m~)/Tf, where Tf is the freeze-out 
temperature. Since, typically, Tf = O(m~/20), this latter factor might suggest that coanni- 
hilation effects would normally be important only for m^ — m~ ~ Tf ~ m~/20. However, a~~ 
is often suppressed by mass and/or phase-space factors in the non-relativistic limit, in which 
case coannihilation processes with larger <7~w, ovy may assume greater relative importance. 
This was indeed found to be the case for coannihilations in the region of MSSM parameter 
space where the LSP is mainly a neutral higgsino that is only slightly lighter than the lighter 

1 However, sneutrinos may be acceptable as dark matter if the MSSM is extended to include additional 
lepton number violating superpotential terms H . 
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chargino and the second-lightest neutral higgsino (TO, [LI 



This is also what we found in a large domain of the MSSM parameter region where the 
LSP is approximately a bino B and the next-to-lightest supersymmetric particle (NLSP) is 
the ?r |T2|]. Moreover, coannihilations with slightly heavier sparticles such as the cr and 
fin are also important. The essential reason is that the non-relativistic threshold S- wave 
contributions to many of the x£r an d £r$*r coannihilation channels are not suppressed by 
fermion mass factors, so that o-±+ ^> <r~~. This is in contrast to higgsino coannihilation 
where ct~h ~ cr~~ above the W ± threshold. 

In a previous paper [12] , we listed many of the important coannihilation channels, re- 
ported on calculations of their cross sections, and emphasized their significance for the 
cosmologically-allowed region of MSSM parameter space. In particular, we highlighted the 
fact that the cosmological upper limit on m~ is increased from ~ 200 GeV, as previously 
estimated |13| , to ~ 600 GeV. This was demonstrated explicitly in the constrained MSSM 
(CMSSM), in which all the soft supersymmetry-breaking mass parameters mo, mi/2 and A&re 
universal at some input supergravity GUT scale. 

In this paper, we amplify and extend this previous discussion by detailing the method 
we have used to calculate coannihilation cross sections and providing simplified analytic 
expressions. These and the numerical results we present serve to explain which coannihilation 
channels are the most important. We also go beyond our previous discussion by discussing 
the dependence of coannihilation effects on such MSSM parameters as A and tan (5. We 
use these results to analyze in more detail not only the impact of coannihilation effects 
on the cosmological upper limit on m~, but also the concomitant bounds on other MSSM 
parameters 0, [H| [16[], such as that on tan/5, which follows in particular from the lower 
limit on the Higgs mass provided by direct searches at LEP. 

The layout of this paper is as follows. In Section 2, we provide a general discussion of relic- 
density calculations, which serves as a framework for our analysis of coannihilation effects. 
Then, in Section 3, we review the standard relic density analysis neglecting neutralino-stau 
coannihilation in the MSSM and CMSSM. In Section 4, we discuss in some detail coannihila- 
tion effects in the CMSSM. Section 5 explores the implications of our coannihilation results 
for the cosmological upper limit on the mass of the LSP and other MSSM parameters, and 
Section 6 draws some conclusions. Useful analytic formulae resulting from our calculations 
are listed in an Appendix. 
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2 Basic Aspects of Relic Density Calculations 



In many cases of interest, the density of relics left over from the early Universe may be de- 
termined relatively simply, once the relevant annihilation cross sections have been calculated 
and used to obtain an annihilation rate. As the Universe expands, a rate or Boltzmann equa- 
tion is solved to determine a freeze-out density, and the relic density subsequently scales with 
the inverse of the comoving volume, and hence with the entropy density. In the case of Dirac 
neutrinos, it is sufficient to calculate an S"-wave cross section to obtain a good approxima- 



tion to the exact result |T7j]. In the MSSM framework discussed here, however, the LSP is a 
neutralino. Since neutralinos are Majorana fermions, the S-wave annihilation cross sections 
into fermion-antifermion pairs are suppressed by the masses of the final-state fermions, and 
it is therefore necessary to compute the P-wave contribution to the cross section || |(| . 
The rate equation for a stable particle with density n is 

^ = - (^rci)(n 2 - n 2 q ) , (1) 

where n eq is the equilibrium number density and (crv re \) is the thermally averaged product of 
the annihilation cross section a and the relative velocity v re \. In the early Universe, we can 
write R/R = {MG N p/?>) 1 / 2 , where p = n 2 g{T)T A /30 is the energy density in radiation and 
g(T) is the number of relativistic degrees of freedom. Conservation of the entropy density 
s = 27r 2 /i(T)T 4 /45 implies that R/R = -f/T - tif/3h where h! = dh/dT. Generally, we 
have h(T) ps g(T). Defining x = T/m and q = n/T 3 h, we can write 

m 



dx 



(^G N g) (h + \mxh>) (av Icl )(q 2 - q 2 J . (2) 



The effect of the h! term was discussed in detail in |18[], and is most important when the 
mass m is between 2 and 10 GeV. Since we only consider neutralinos that are significantly 
more massive, we neglect it below. 

For x ^> 1, neutralinos are in thermal and chemical equilibrium and q = q eq ~ constant, 
since n eq ~ T 3 . When x = 0(1), q — q eq ~ e~ x l x jx z l 2 until freeze-out, after which q is 
again approximately constant. For annihilations governed by weak-strength interactions, 
freeze-out occurs when x ~ 1/20. The final relic density is determined by integrating (|I]) 
down to x = 0, and is given by 

p x = mq(0)h(0)T* (3) 

More generally, when coannihilations are important, there are several particle species i, with 
different masses, and each with its own number density rij and equilibrium number density 
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n. 



cq 



i j. In this case JTIJ, the rate equation ([I]) still applies, provided n is interpreted as the 



total number density, 



n 



E 



the total equilibrium number density, 



n, 



cq 



E 



n, 



eq,i ; 



(4) 



(5) 



and the effective annihilation cross section as 

^eq,i^eq,j 



(a eS v ie i) = 



,, "cq 



(6) 



In eq. (0), m is now understood as the mass of the lightest particle under consideration. 
For T < rrii, the equilibrium number density of each species is given by JT^, [19|] 



n 



eq,i 



9i 



d 3 p 
(2^)3 



-E/T 



9l m\T 



K 2 (m l /T) 



9i 



2tt 2 



2tt 



J exp(-m;/T) ^1 + 



15T 



sm,; 



(7) 



where ^ is a spin degeneracy factor and K 2 (x) is a modified Bessel function. We have made 
the approximation of Boltzmann statistics for the annihilating particles, which is excellent 
in practice. 

We now wish to compute (cr^frei) for the process 1 + 2— >3 + 4inan efficient manner. 
Suppose we have determined the squared transition matrix element |T| 2 (summed over final 
spins and averaged over initial spins) and expressed it as a function of the Mandelstam 
variables s, t, u. We now wish to express |T| 2 in terms of s and the scattering angle 8 CM in 
the center-of-mass frame. We have 

(m 2 — m^mj — m 2 ) 



t — u 



+ 4p 1 (s)p 3 (s) cos d CM , 



where Pi{s) is the magnitude of the 3-momentum of particle i in the CM frame, given by 

1/2 



Pl[S) =P2{S) 



P3{S) = Pl{S) 



s ml + m\ (m 2 — m 2 ) 2 
4 2 + As 



s m\ + m\ (m 2 — m|) 2 
4 2 + Is 



1/2 



(9) 
(10) 
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We can then use s + t + u = m\ + m 2 + m| + m 2 to write 



t = \ [m 2 + m 2 + m 2 + ml — s + (t — u)] 
u — 2 [ m i + m 2 + m 3 + m 4 — s — (t — n)] 

Using (|8|)-(ll2D, we can write |T| 2 as a function of s and cos# CM . 
We now definef] 



(11) 
(12) 



w(s) 



d 3 p 3 d 3 p 4 



4 J (2tt) 3 £ 3 (2it) 3 E± 

1 p 3 ( S ^ r+l 



(27r) i 6 i (p 1 +p 2 -p 3 -p 4 ) \T\ 



32tt gi/a 7_ 



(icos^M T 



(13) 



In terms of w(s), the total annihilation cross section oi 2 (s) is given by 012(5) — w(s)/ s l ^ 2 p\{s). 
Our ty(s) is also the same as w(s) in ]Ll| [|, [12]], which is written as W/4 in |20] . 

So far all this is exact. To reproduce the usual partial wave expansion, we expand |T| 2 
in powers of px(s)/mx. The odd powers vanish upon integration over 8 CM , while the zeroth 
and second order terms correspond to the usual S and P waves, respectively. We see from 
(H) that each factor of p\(s) is accompanied by a factor of cos# CM , so we have 



+1 



dcos8 CM |T 



T 



2 

cos 6»cm ■ 



'+1/V3 



+ iri 



cos^cm — » -l/v 7 ^ 



(14) 



We can therefore evaluate the S and P wave contributions to w(s) simply by evaluating |T| 2 
at two different values of cos# CM ; no integrations are required. 



The proper procedure for thermal averaging has been discussed in [|18| , [IQf] for the case of 
mi = m 2 , and in [|, |2(J for the case of mi 7^ m 2 , with the result 

1 



(Cl2^rel) 



2m 2 l mlTK 2 (m 1 /T)K 2 (m 2 /T) . '(mi+mj) 



dsXi(x/i/T)pi(s)w(s) . (15) 



Using the asymptotic expansion K n (x) = (2x/7r)- 1 / 2 e^'[l + (4n 2 -l)/(8x) + ...] of the Bessel 



functions, changing the integration variable from s to y 
expanding in powers of T, we find 

w(s ) 3 (mi + m 2 ) \w(s ) 



,1/2 



mi — m 2 )/T, and then 



(0-12^rel) 



m 1 m 2 
1 

m\m 2 



2m\m 2 



mim 2 



2w'{s ) 



T + 0(T 2 



3(mi + m 2 )T 



2m\m 2 



w(s)\ 



s~ t(mi+m 2 ) 2 +3(mi+m2)T 



+ (9(T 2 



«12 + &1 2 X + C(X 2 



(16) 



2 An extra factor of i must be included in the case of identical final state particles. 
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where x = T/m\ (assuming ni\ < 777.2), and, in the first line, So = (mi + m2) 2 . We extract a±2 
and 612 from the transition amplitudes listed in the Appendix by performing the substitutions 
(§)-([L6|) for each final state. We set x = to get a i2 , and then compute 612 by setting x 
to a numerical value small enough to render the 0(x 2 ) terms negligible. We compute a e g 
and b e ff by performing the sum over initial states as in eq. (|B|). We then integrate the rate 
equation @ numerically to obtain the relic LSP abundance. To a fair approximation, the 
relic density can simply be written as || 10 



ah > „ "■•<**- (17) 

g/ M p i(a eS + b eS x f /2)x f 

where the freeze-out temperature 2} ~ m~/20, and g f is the number of relativistic degrees 
of freedom at 7}. Note that this implies that the ratio of relic densities computed with and 
without coannihilations is, roughly, 

where a = a + bx/2 and sub- and superscripts denote quantities computed ignoring coanni- 
hilations. The ratio x®/x f ~ l+x® \n(g e ffa e e / giao) , where g e s = Y,i gi( m i/ m i) 3 ^ 2 e~^ mi ~ mi ^ T ■ 
For the case of three degenerate slepton NLSPs, g c s = Ya 9% = 8 an d x ?/ x f ~ 1-2- 

The non-relativistic expansion (|I^) is known to be inaccurate near s-channel poles and 



final-state thresholds [0, [19], where the cross section can depend strongly on the initial 
LSP momenta, and where LSPs in the Boltzmann tail can access resonances and final states 
forbidden to zero-momentum LSPs. In the CMSSM, well-studied examples are neutralino 
annihilation on the Z and light Higgs poles, where a detailed treatment of the thermal aver- 



aging [Ty, [19j is necessary to compute accurately the neutralino relic abundance. Examples 
of final-state thresholds include it and hZ. 

However, none of these effects are significant for the analysis of this paper. The regions 
with my, < mz/2 and m% ^ ^/i/2 have now been excluded by LEP chargino searches. 
Annihilation on the H and A poles can be important for tan (3 > 40 [21] , where the heavy 
Higgs masses can be close to 2m^ for very large ranges of my 2 (since the heavy Higgs and 
bino masses scale similarly with myz). However, we do not consider in this paper such large 
values of tan/3, where stau mixing may be important as well. In the CMSSM, the only 
non- negligible threshhold occurs at m^ ~ m t , which is visible as a kink near my 2 = 400 GeV 
in Figs.[l|4J. However, as can be seen there, the effect of the top threshhold (and hence also 
that of sub-threshhold neutralino annihiliation into tops) on the total annihilation rate is 
tiny, because the contribution from the it final state is suppressed by the large stop masses. 
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Likewise, the H and A masses are much too large in the CMSSM for heavy Higgs final 
states to be of relevance, except perhaps again for the very large values of tan (3 > 40 (not 
considered here), where the heavy Higgs masses are smaller. Finally, other thresholds such 
as hZ are suppressed by the smallness of the Higgsino component of the LSP, which is a 
very pure bino in the CMSSM. 

As for such effects in coannihilations (via sleptons in our case), we find that coannihila- 
tions are important in the CMSSM for fairly large values of my 2 > 300 GeV. The relevant 
CMSSM regions are thus above the light thresholds (e.g., £^£^ — > W + W~) and still well be- 
low the heavy threshholds (e.g., £^£^ — > 'jH) for slepton annihilation and slepton-neutralino 
coannihilation, since only staus with masses rrif R ~ m% affect the neutralino relic density. 
The process £^ £^ — > tt is numerically irrelevant (see Figs.|l]-[|]). Therefore the partial- wave 



expansion (14) is a valid approximation in our analysis. 



3 The MSSM Without Coannihilations 

In the MSSM, the identity of the LSP is determined by the following tree-level parameters: 
the supersymmetry-breaking gaugino mass m ly / 2 (assuming gaugino mass universality at 
the GUT scale), the supersymmetric Higgs mixing mass //, and the ratio of Higgs vacuum 
expectation values (vev's), tan/3. The annihilation cross section and hence the relic density 
depend on the identity of the LSP. For mua the LSP is mostly a higgsino state, 

whilst in the opposite limit, /i ^> m^, the LSP is mostly a bino B. The negative results 
of previous SUSY particle searches at LEP have been able to impose strong constraints on 
the MSSM parameter space 0, [15|, [l|J. Roughly speaking, one may conclude that m\/2 and 



H > 80 GeV. 

In the region of the MSSM parameter space where the LSP is mainly a B, annihilations 
to fermion-antifermion pairs proceed mainly through sfermion exchange. As we discuss in 
more detail below, this process is P-wave suppressed, which implies a great reduction in the 
annihilation rate at the low temperatures at which the P's freeze out. Over much of this 
region, it is possible to obtain a cosmologically significant relic density, as we review below. 
By contrast, in the region where the LSP is mainly a higgsino, annihilations to W ± pairs 
are dominant above threshold. This process is not P-wave suppressed, and, as a result, the 
relic density is suppressed in this higgsino region, unless one considers either very heavy 
higgsinos (mg > 500 GeV), or higgsinos with masses below the W threshold. Further, in 
much of the higgsino region of parameter space, either the second-lightest neutralino or the 
higgsino-like chargino is very close in mass to the LSP. In this case, it was shown |TT[ 
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that coannihilations were important in determining the relic density of LSPs. Indeed, when 
combined with the current LEP limits, there remains very little (if any) parameter space 
remaining where a light higgsino LSP could contribute a sufficient relic density (Qh 2 > 0.1) 
to be of interest as the dominant component of astrophysical dark matter At very 

low tan/? (1.2-1.6), some solutions with Qh 2 > 0.1 may exist just below the W threshold 
|2(| . However, the LEP Higgs bounds (though dependent on quantities such as mi) make it 



extremely difficult to achieve tan/3 this low, as discussed, for example, in [[LJ, [T5|, [lq| . 

It is also of interest to consider a constrained version of the MSSM (CMSSM), in which 
all the soft supersymmetry-breaking scalar masses, mo, are unified at the GUT scale. In 
this case, the conditions which determine electroweak symmetry breaking also fix and 
the pseudoscalar MSSM Higgs mass at the weak scale. For all choices of m\/2 and m 
consistent with LEP mass bounds, the lightest neutralino is predicted to be a B, modulo thin 
fringe strips of parameter space close to where the electroweak symmetry is not dynamically 
brokenf]. Although B's typically have an interesting relic density, this is no longer true if the 
B mass happens to lie near m^/2 or m/j/2, in which case there are large contributions to 
the annihilation through direct S-channel resonance exchange. However, since LEP limits 
on the chargino mass can be translated into bounds on mi/2, these resonant cases are now 
all but excluded for small values of tan /3, as we discuss further in Section 4. 

Since we are mainly interested the B LSP candidate, we focus our discussion on this case, 
in order to be more specific. The thermally-averaged cross section for BB — > ff takes the 
generic form 



(av) - 



1287T 




1/2 



(Yl + Y 2 R?['^)+ V? + Y£) { ^4) (l + 0(m 2 /m 2 s ))x 




(19) 

where Yl{r) is the hypercharge of Jl(R), Ay = m| + m~ — m 2 , and we have shown only the 
leading P-wave contribution proportional to x = T/mg. As advertised, the S"-wave piece 
is proportional to the fermion mass-squared and hence is negligible, except perhaps for the 
top quark, and this has the net effect of reducing the neutralino annihilation cross-section 
by 0(xf). The general form of ( |i~9|) leads to an upper bound on the possible mass of the 
LSP within the MSSM, due to the cosmological relic density |13j. Specifically, in the case 
of a bino LSP, the upper limit on comes about as follows. The assumption that the B 
is the LSP requires, in particular, that mg < mj. In order to minimize the relic density, 
we must maximize the cross section, which is done by setting mj = m~. The cross section 

3 In fact, unification of the gaugino masses at the GUT scale by itself implies a S-like LSP over the bulk 
of the parameter space [||| . 
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is then approximately inversely proportional to m~. The cosmological upper limit on Q^h 2 
translates into a lower limit on (av) which then, in turn, yields an upper limit to mg. In 
the MSSM, this limit is mg < 300 GeV, when all sfermion masses are taken to be equal at 
the weak scale, though the limit can be weakened when sfermion mixing [23| or CP-violating 
phases are included [E3 . 



In the CMSSM, the argument is somewhat similar, although mg and the sfermion masses 
are now no longer entirely independent, because it is assumed in the CMSSM that there is 
a common scalar mass m at the GUT scale. For a given value of the common gaugino 
mass my? at the GUT scale, the relic B density falls with decreasing mo, since 



m f = 

m o + C/ m i/2 + C( m l); where Cf is a positive numerical coefficient that is calculable via 
the renormalizat ion-group evolution of the sfermion masses. Therefore, the cosmological 
upper limit on Q^h 2 translates at fixed into an upper limit on m . At low to moderate 
tan/3, this upper limit is typically not larger than mo ~ 170 GeV, unless one is sitting on a 
direct-channel pole, i.e., when mg ~ m^/2 or in which case s-channel annihilation is 

dominant and there is no upper limit on m . However, we are interested in an upper bound 
on mg, and hence in masses far from the light Higgs and Z poles. We recall that scales 
with my 2, and it transpires for vnyi > 400 GeV that mg exceeds the mass of the lightest 
sfermion, which is typically the tr, for m small enough to satisfy the cosmological bound 
as traditionally computed (i.e., neglecting coannihilations) P5| . Thus, the LSP is no longer 
a neutralino for such large values of m^, and hence an upper bound mg < 200 GeV P5| 
can be established. ^ 

In Jl5| , it was shown that the LEP constraints on the mass of the supersymmetric Higgs 
boson, when combined with the above cosmological upper limit on the LSP mass (or m^), 
leads to an interesting bound on tan/3. The argument is as follows. At fixed tan/3, due 



to the radiative corrections to the Higgs mass in the MSSM p9 |, it is always possible to 

satisfy a given experimental constraint on the lightest MSSM Higgs mass by going to large 

values of either or tuq. At large mo, however, the cosmological bound forces one to low 

so that the LSP is higgsino-like and annihilations are not suppressed by large sfermion 

masses. In the CMSSM, this is not possible, since /i is fixed by the condition of electroweak 

symmetry breaking. Even when universality among the soft Higgs masses is not assumed, if 

m is large, obtaining a reasonable relic density requires to be low enough to ensure that 

the LSP is a higgsino. However, unless m 1( / 2 is very large, lowering this far results in a 

lower Higgs mass. At a given value of tan /3, the Higgs mass bound can be translated into 

4 This upper bound can be strengthened by requiring that the global minimum of the effective potential 
of the MSSM conserve electric charge and color Eg, . 
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a lower bound on mu2. If this lower bound is greater than the cosmological upper bound 
on mi/2 discussed above, the corresponding value of tan/5 considered is excluded. Since the 
lower bound on mx/2 due to the Higgs mass bound is very dependent on tan /3, we can derive 
a lower bound on tan (3 by combining the LEP and cosmological bounds |T3], . 

4 Coannihilations in the CMSSM 

As discussed earlier, if the masses of the next-to-lightest sparticles (NLSPs) are close to the 
LSP mass: AM = 0(xj)M, the number densities of the NLSPs have only slight Boltzmann 
suppressions with respect to the LSP number density when the LSP freezes out of chemical 
equilibrium with the thermal bathf]. Moreover, it is well known that, in such circum- 
stances, coannihilations of NLSPs with the LSP, along with NLSP-NLSP annihilations, may 



play an important role in keeping the LSPs in chemical equilibrium with the bath ||10|| , and 
the number density of neutralinos can be significantly reduced by such coannihilations. These 
processes can be particularly important when the LSP annihilation rate itself is suppressed, 
as is the case for neutralinos, as discussed above. The case of heavy higgsinos is a well 
studied example f50fl . Analogously to that case, the B relic density can be reduced through 



coannihilation with slightly heavier r R 's or other sleptons, as we now discuss in detail. 

In the CMSSM, when attains the upper bound discussed in the previous section, the 
B is degenerate in mass with the f R , and quite close in mass to the and Jjlr. In [|12| , 
we showed that the effects of coannihilations between the neutralino LSP and the tr (also 
including the jiR and Sr) can have have a dramatic effect on the derived upper bound on 
mi/2 and the mass of the LSP. Such coannihilation effects thereby also affect the derived 
bound on tan j3. 

To compute the effective annihilation cross sections for light sparticles in the CMSSM, 
we allow the index % in (f|) to run over f R , f*, e R , e* , /i R and /I*, as well as %■ Many of the 
resulting 49 in (H) are related, so we can write 

Ocff = o xx r x r x + 4 a XT r x r T + 8 o xe r x r e + 2 (o TT + cr TT *)r T r T + 8 (a re + a Te *)r T r e + 

4 (cr ee + a ee *)r e r e + 4 (a eM + a e ^*)r e r e , (20) 

where = n eq ^/n eci , we have taken the e R and J1 R to be degenerate in mass, and we have 

neglected the electron and muon masses. We list in Table [I] the sets of initial and final 

5 We recall that 2-2 scatterings with particles in the thermal bath keep the NLSPs in chemical equilibrium 
with each other and with the LSP, down to temperatures well below the temperature at which the comoving 
LSP number density freezes out. 
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states for which we compute the annihilation cross sections. We list in the Appendix the 
transition amplitudes for the scattering processes, which are sufficient to compute the a 
and b coefficients, following the discussion in Section 2. We have verified that the r mass 
is numerically irrelevant in our analysis, and the formulae we present are simplified to the 
m T —>■ limit, although our numerical results do include a non-zero r mass, including extra 
diagrams which are present only when m T ^ 0, but which are not listed in the Appendix. 
We ignore the effect of f mixing on the cross-sections, which may be important at large 
tan f3. The final states involving heavy Higgses are kinematically unavailable in the regions 
of CMSSM parameter space relevant to our analysis. However, for completeness, their tran- 
sition amplitudes are also tabulated in the Appendix. As %f coannihilation is important at 
large my2, where the B purity is very high, we compute the three xj amplitudes in the B 
limit, where t-channel neutralino exchange is suppressed. We also have not included three 
body final state processes, such as s-channel Higgs exchange to a Higgs + gauge boson pair, 
via a two Higgs/two gauge boson vertex. In addition to being phase-space supressed, the 
magnitude of these contrubtions to the total cross section is down by a factor of a relative 
to the two body final states we have included. 



Table 1: Initial and Final States for Coannihilation: {i,j = T,e,fi} 



Initial State 


Final States 


7i pi* 
R C R 


77, ZZ, jZ, W+W-, Zh,jh,hh, ff, 




ZH, 7 H, ZA, W + H- hH, hA, HH, HA, AA, H+H 


C R l R 




^R ^R 5 * 7^ 3 






i^.tZ.th 



We display in Fig. [I] numerical values of the contributions to a = a + bx/2 (see fllTD), 
for the representative values mo = 120 GeV, Aq = 0,x = 1/23, tan/? = 3, and fi > 0, as 
a function of my 2- F° r comparison, the total cross section for XX annihilation to all final 
states is shown as a thick dotted line. Due to the P-wave suppression of the cross section 
for xx annihilation to fermion pairs, the XX cross section tends to be an order of magnitude 
smaller than the others, which is why coannihilation effects are so important. In practice, 
we find that the dominant contributions to a c g generally come from annihilations of £^£^ 
to gauge bosons, £^ to lepton pairs, and £^ % to £} + gauge boson. Due to the momentum 
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200 300 400 500 600 700 800 
m 1/2 

Figure 1: The separate contributions to the cross sections a = a + ^bx forx = T/m~ = 1/23 
and m = 120 GeV ; as functions of my 2 : a)ff*, b)fx, and c) other interactions. For 
comparison, the thick dotted line is the XX cross section. 



12 



CD 

o 



10- 




10- 



'10- 



T 

IX 
it- 



<b 



200 300 400 500 600 700 800 
m, 



L l/2 



tan/3 = 3 



G 



TOT 



In I I ''H'hrvi. >,i.,, i i h 

200 300 400 500 600 700 800 
m 1/2 



10- 



> 
■-. 
a 



40- 



tari(S = 3 



G 



TOT 



In mm rri-i-ij I 'i'H'hrvt.tii.,,1 i h- 

200 300 400 500 600 700 800 
m 1/2 



Figure 2: j4s m Fzg. Q but for the choice m = 200 GeF. 
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m 1/2 

Figure 3: As in Fig. [7], but for the choice tan/3 = 10. 
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200 300 400 500 600 700 800 

m i/2 

Figure 4: As in Fig. but for the choice tan/3 = 10, and A = —3my2- 
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dependence of the ff*Z and ff*7 couplings, the cross sections for ff* annihilation into 
Zh, 7/1 and // are also P-wave suppressed. The Zh and 7/1 cross sections are thus off the 
bottom of the figure, and // is significant only by dint of the large number of final states. 
The top threshold is visible as a small bump on the // line between my 2 of 300 and 400 GeV. 
The curve for ff* — > W + W~ in the first panel is obscured by the thick dotted line. As 
can be seen by comparing Figs, [l] and |2|, there is a dependence of the cross-sections on 
m , and the variation is largest at low ray 21 where the slepton masses are most sensitive to 
m . The sharp rise near my 2 = 200GeV in Fig. ||]a is due to heavy Higgs poles. Figures for 
tan (3 = 10 are very similar to those for tan (3 = 3, as can be seen in Fig ||. There is hardly 
any dependence on A for tan/3 = 3, whilst, for tan/3 = 10, a mild dependence can be seen 
in Fig |], where we have taken A = —3my 2 . 

The contributions of the various annihilation channels to a e g are weighted by the rel- 
ative abundances of the tr, Sr, Jar and \- F° r sleptons degenerate with the %■> slepton 
annihilation and slepton-neutralino coannihilation clearly dominate the contributions to 
cr e ff in (H), and the final neutralino relic density is greatly reduced. As the sleptons be- 
come heavier than the neutralinos, their number densities are exponentially suppressed, 
and when the mass differences are no longer small, the slepton contributions to a c g are 
negligible. Fig. |5| shows the sizes of the separate contributions to a e g from neutralino an- 
nihilation, neutralino-slepton coannihilation and slepton-slepton annihilation and coannihi- 
lation, as functions of the mass difference between the f and x- m Fig- |§ a > we have fixed 
my 2 = 500 GeV, tan (3 = 3, A = 0, fi > 0, and computed a e g for varying m , which amounts 
to varying the mass difference AM. In this case, all the sleptons are closely degenerate with 
each other and hence contribute equally to a e s. The thin solid line is the a which one would 
compute if one ignored the slepton states, i.e., a xx + b xx x/2. Note that, in the case of close 
degeneracy between the \ and tr, it is in fact slepton annihilation by itself which dominates 
(j e fj. Since this contribution is suppressed by two powers of n^^, it drops rapidly with 
AM, and neutralino-slepton coannihilation takes over at AM > 0.07. This contribution in 
turn falls with one power of n eqi f R , and neutralino annihilation re-emerges as the dominant 
reaction for AM > 0.17. When AM > 0.25, the two solid lines and dotted line all merge, 
and coannihilation can be neglected. Figures for other my 2 and tan f3 are shown in Fig. [5|b-d: 
the underlying physics is similar, and they are similar in shape to Fig. |^a, but with differ- 
ent normalizations. In Fig. [5|b, we take tan/? = 10, in Fig. |5]c,d we take tan/? = 3, with 
m i/2 — 300 and 1000 respectively. The values of m x shown on the figures is determined by 
the choice of my 2 . 
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Figure 5: The separate contributions to the cross sections a e g for x = T /m~ = 1/23, as 
functions of AM = (mf R —m~)/m~, with a) (mi/2,tan/3) = (500 GeV, 3), b) =(500 GeV, 10), 
c) =(300 GeV, 3), and d) =(1000 GeV, 3). 
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5 Implications of Coannihilations for the Upper Limit 
on the LSP Mass 



We now explore the consequences of coannihilation for cosmological bounds in the CMSSM. 
In Fig. we display the cosmologically and experimentally permitted regions of the {my 2, mo) 
plane. We have chosen the representative points tan/3 = 3 and 10, and present results for 
both fi < and \x > 0. The light shaded regions correspond to 0.1 < ^h 2 < 0.3. The dark 
shaded regions have m~ R < m~ and are excluded by the very stringent bounds on charged 
dark matter 0. The light dashed contours indicate the corresponding regions in ^^h 2 if 
one ignores the effect of coannihilations. Neglecting coannihilations, one would find an upper 
bound of ~ 450 GeV on m^, corresponding to an upper bound of roughly 200 GeV on Wg. 
The effect of coannihilations is to create an allowed band about 25-50 GeV wide in mo for 
myi < 1400 GeV, which tracks above the m fR = m~ contour. Along the line = m~, we 
find R ~ 10, as shown numerically in Fig. [5] and fli~8|) . As mo increases, AM increases and 
the slepton contribution to cr cff falls, as in Fig. |, and the relic density rises abruptly. 

We also display bounds from LEP particle searches [51] in Fig. |6|. The chargino mass 
bound from LEP essentially saturates the kinematic limit of ~ 95 GeV, modulo a small 



gap which occurs when the chargino is just slightly heavier than the sneutrino |3T[. The 



chargino 95 GeV isomass contour is displayed as a dark dashed contour, and it cuts off most 
of the annihilation pole zone, which appears as a chimney in the cosmologically preferred 
regions at low my 2 for the values of tan f3 considered. Representations of slepton limits from 
searches |3jJ for acoplanar lepton pairs at LEP 183 appear as light dotted contours. The 
most severe experimental constraint at low tan /3 comes from LEP Higgs searches. At low 
tan p, the tree level light Higgs mass m^ mz\ cos 2/3\ lies well below the experimental bound 



> 95 GeV. Radiative corrections to rrih are known to be large |2{| and increase logarithmically 

with the stop masses, and hence with my%. Thus, for sufficiently large my 2, the Higgs bound 

may be satisfied. However, if the minimum my 2 exceeds its cosmological upper bound for 

a given tan (3, this value of tan f3 is excluded. In the absence of coannihilations, the lower 

bound on tan/? was computed |L6|] to be 2.0 for /x > and 1.65 for \i < 0, using mass 

limits from LEP 183. Higgs isomass contours of 95,100,105, and 110 GeV are displayed as 

dot-dashed lines in Fig. |6]. 

When we include coannihilations, the cosmological upper bound on my 2 weakens con- 

6 To be more precise: here and subsequently, we have included consistently the effects of tr — ft mixing on 
the mass of the lighter f (which is mainly tr), both in delineating the cosmological exclusion domain and in 
the kinematics of coannihilation. However, we have not included mixing angle effects in the (co) annihilation 
amplitudes, since these are small for the values of tan/3 studied in this paper. 
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Figure 6: The light-shaded area is the cosmologically preferred region with 0.1 < ^^h 2 < 0.3. 
The light dashed lines show the location of the cosmologically preferred region if one ignores 
coannihilations with the light sleptons. In the dark shaded regions in the bottom right of each 
panel, the LSP is the tr, leading to an unacceptable abundance of charged dark matter. Also 
shown are the isomass contours m x ± =95 GeV and = 95, 100, 105, 110 GeV, as well as 
an indication of the slepton bound from LEP In the area below the solid contour, the 

scalar potential contains charge and/ or colour breaking minima. 
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siderably, as we have already mentioned. In Fig. [7], we extend the my 2 coverage to larger 
scales, to show the cross-over point between the regions with VL~h 2 < 0.3 and m~ < m~. 
The two constraints together require my 2 ^ 1450, corresponding to an upper bound on the 
neutralino mass of m~ < 600 GeV. We take 93 GeV as our experimental lower limit on the 
Higgs mass, to allow for a 2 GeV uncertainty in the extraction of the radiative corrections 
to the Higgs mass [29|. The combined lower bound on tan/? from the Higgs search and 
cosmology, including the recent LEP 189 limits, indicates that tan/? > 2.2 for \i < and 
> 1.8 for fj, > Q. 

Figs. H and |7] were generated for the particular choice A = —my 2 , for reasons described 
below. Since the stau mass depends weakly on A , particularly at larger values of tan/?, 
the position of the m,f R = m~ contour shifts somewhat with A . In Fig. [| we show figures 
corresponding to Figs. |6]a and ||c for A = —3my 2 . For tan (3 — 3, the width of the cosmolog- 
ically allowed band is the same as in Fig. |]a, but the band becomes smaller for tan (3 = 10. 
The effect for positive Aq is qualitatively similar. Likewise, varying Aq has a tiny effect on 
the cosmological upper bound on my 2 for tan/? = 3, but can produce a small (~ 100 GeV) 
reduction in the upper bound for tan/? = 10. 

Lastly, we show as a dark solid contour the constraint coming from the requirement that 
the global minimum of the scalar potential not break charge and colour (CCB) P6] , PS] , f2"T| . 
The areas below the light solid line in Figs. |G] and [7] contain charge and/or colour breaking 
minima, while the regions above the lines are free of such minima^, modulo a thin (~ 10 GeV 
wide) strip on top of the solid contour, where local (but not global) CCB minima exist P71 . 
We have chosen A = —my 2, where the CCB bounds are weakest ||26|| . The bounds are not 
strongly dependent on the sign of /z and are strongest for low tan /?, where the top Yukawa 
coupling is largest. In the absence of coannihilations, it is clear from Fig. ^ that a fair 
fraction of the cosmologically allowed regions contain charge and colour breaking minima. 
In the presence of coannihilations, the cosmologically allowed region is extended, and as the 
CCB eventually falls with my 2-, CCB free regions with small relic densities may occur in the 
proboscis region as well, particularly at the larger values of tan/? considered. In Fig. [7|a, 
we see that, for tan/3 = 2, the entire cosmologically allowed region of parameter space also 
contains CCB minima, even for the conservative choice A = —my 2 . For larger tan/?, the 
CCB contour bends over and, for Aq = —my 2, exposes a small strip of the shaded region 
above the m fR = m~ line. For comparison, we also show the corresponding CCB contours for 

7 The LEP bounds on the CMSSM parameter space will be considered in more detail in |}2j , as well as a 
comparison with the supersymmetric reach of the Fermilab Tevatron collider. 

8 The difference between the CCB curves of Fig. || and |27j is due to our choice here of m t — 175 GeV, 
rather than 170 GeV. 
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Figure 7: The same as Fig. for fi < 0, but extended to larger m^. The dashed contours 
in panels b) and c) correspond to the solid CCB line for A = 0,—2my2, and the dotted 
contour for m t = 170 GeV. 
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Aq = — 2m 1 / 2 and 0. For these choices of Aq, there is again no region which satisfies both the 
CCB and cosmological constraints. The CCB bounds are also sensitive to the top mass. In 
Figs. ^ and |7| we have taken m t = 175 GeV; taking m t = 170 GeV produces the dotted CCB 
curves in Fig.^a-c, which expose a much larger piece of the cosmologically allowed band. We 
recall, however, that a low value of m t also reduces the radiative corrections to the Higgs 
mass, and therefore makes the Higgs constraint more severe. However, this does not cause 
a conflict with the cosmological constraints, even for tan/3 as low as 3. 

There is one caveat which one must bear in mind when interpreting the CCB bounds: 
the tunneling rate from the charge and colour conserving minimum to the CCB minimum 
is very slow, to such an extent that the conserving minimum is essentially stable over the 
current lifetime of the universe. Therefore the presence of CCB minima may present more of 
a cosmological problem of how to populate the physical minimum in preference to the CCB 
minimum than a constraint on the particle physics |33] . 



We expect qualitatively similar effects on the corresponding bounds in the MSSM, that is, 
when we drop the condition of universality of scalar mass at the GUT scale. Though there is 
an upper limit of 0(a. few TeV) for Higgsinos in the MSSM, the limits due to coannihilation 
that we have been concerned with apply only to the gaugino limit. In this case, one must 
in general take all the squarks and sleptons degenerate with the neutralino and compute 



22 



the annihilation and coannihilation cross sections for all possible combinations of sfermions. 
However, if the rates are the same as for the sleptons, the effect is about a 10-15% decrease 
in (Q~h 2 )°/R, leading to a similar bound on m~ as in the CMSSM. In practice, of course 
the result will depend on which state is the NLSP. 



6 Conclusions and Open Issues 

We have documented in this paper the importance of coannihilation effects on the relic den- 
sity of B dark matter. We have presented a general formalism for including their effects in 
relic density calculations, and have provided (in the Appendix) analytic approximate formu- 
lae for many coannihilation processes involving tr, cr and fiR sparticles. We have given in 
the Figures numerical examples that exhibit their relevance for various values of the CMSSM 
parameters mi/2,fn , A and tan/3. We have explained why coannihilation effects are so im- 
portant, principally because of the P-wave suppression of non-relativistic \X annihilation. 
One immediate physical consequence of these coannihilation effects is to relax signifi- 



cantly the previous cosmological upper bound [13| on the LSP mass, from > 200 GeV to 
> 600 GeV. This relaxed upper limit could have significant implications for strategies to 
search for supersymmetric cold dark matter, which have yet to be explored systematically. 
Coannihilations also allow the LSP mass to approach the boundary of the CMSSM param- 
eter space within which LHC searches for supersymmetry are expected to be sensitive |§. 
At first glance, it seems likely that the LHC should still be able to cover all of the CMSSM 
parameter space that is consistent with supersymmetric cold dark matter, but this point 
merits further consideration. 

Another point to be studied in more detail is the interplay between coannihilation effects 
and LEP lower limits on sparticle and higgs masses. Previously, limits on m~ and tan/5 have 



been derived from a combination of previous LEP data sets and cosmology |L4|, [H| [T(| , but 
neglecting coannihilation effects. These limits should now be revisited in the light of more 
recent LEP data |3T|, as well as coannihilations. We have commented on these questions in 
this paper, but a complete study lies beyond the scope of the present analysis. 

We hope that this paper has not only documented the importance of coannihilation 
effects, but also provided the reader with the tools needed to join in the exploration of their 
implications for the interesting open issues raised in the previous paragraph. 
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Appendix 

This section contains simplified formulae for the ff*,JX, ff,f£ and ft annihilation ampli- 
tudes, in the m f — > limit. Expressions for the ££*,££ and £\ amplitudes can be obtained 
by taking r — > I in the ff*, fx, ff formulae below, and fie and ejl* by taking r — > e, £ — > fi 
in the fi and fi* expressions. The a and b coefficients can be simply numerically extracted 
from the amplitudes, as described in section 2 of the text. 

ff* — ► w+w- 

I. s-channel H annihilation 

II. s-channel h annihilation 

V. s-channel Z annihilation 

VI. s-channel 7 annihilation 



fi = (— g2in w cos(P — a)) (g 2 m z sin 2 9 W cos(a + P)/ cos 9 W ) 

fi = (— g 2 mw sin(P — a))(—g2m z sm 2 9 w sm(a + (3)/ cos9 w ) 

h = (-92 sin 2 9 w / cos 9 w )(g 2 cos 9 W ) 

h = e 2 

71x71 = {I2m 4 w - 4m 2 w s + s 2 )/(4m 4 w {m 2 H - s) 2 ) 

T n xT u = (12m 4 w -4m 2 w s + s 2 )/(4m 4 w (m 2 h -s) 2 ) 

TyxTv = (128m 2 R m 4 v m 4 z s — 32m 2 R m 2 v m 4 z s 2 — 32m 4 v m 4 z s 2 + 8m 2 v m 4 z s i + 12m^r?f|t 2 — 

ey a c\ /IOO A A 'y A A 'y A A o 

4m w m z st + m z s t — 24m w m z tu + 8m w m z stu — 2m z s tu + 12m w m z u — 
4m 2 v m z su 2 + m z s 2 u 2 ) / \4m 4 v m 4 z (m 2 z — s) 2 ) 
T YI xT yi = (128m 2 fR m 4 v s - 32m\ R m 2 w s 2 - 32m 4 v s 2 + 8m 2 w s 3 + \2m 4 w t 2 - 4m^st 2 + 
s 2 t 2 - 2Am\ytu + Smlystu - 2s 2 tu + 12m^w 2 - 4m 2 y sw 2 + s 2 u 2 )/(4m 4 v s 2 ) 
TixTu = (12m 4 v -4m 2 w s + s 2 )/(4m 4 w (m 2 H -s)(m 2 h -s)) 

%xT y = (— I2m 4 v m 2 z t + m 2 z s 2 t + 12m 4 v m 2 z u — m 2 z s 2 u)/(4m 4 v m 2 z (m 2 H — s){m 2 z — s)) 
TiixTy = (— Y2m w m z t + m z s t + 12m w m z u — m z s u)/(4m w m z (m h — s){m z — s)) 
TjxTyj = -{-\2m 4 w t + s 2 t + \2m 4 w u - s 2 u) / l (4m 4 v (m 2 H - s)s) 
T u xTyi = -{-\2m 4 w t + s 2 t + \2m 4 w u - s 2 «)/(4m^(m^ - s)s) 
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TyxTyi = [—\28m\ R m\ v m 2 z s + 32mf R m 2 v m 2 z s 2 + 32m 4 v m 2 z s 2 — 8m 2 v m 2 z s 3 — 
12m\ym 2 z t 2 + Am 2 w m 2 z si 2 — m 2 z s 2 t 2 + 2Am\ v m 2 z tu — 8m 2 v m 2 z stu + 
2m 2 z s 2 tu — 12m 4 v m 2 z u 2 + 4m 2 1/ m|s-u 2 — m 2 z s 2 u 2 ) / (Am\ym 2 z (m 2 z — s)s) 
\T\ 2 = f 2 T 1 xT 1 + f 2 T n xT u + f 2 T v xT Y + / 6 2 T V ixT V i + 2f 1 f 2 T l xT ll + 2/ 1 / 5 ^xT v + 

2/i/ 6 ^xT VI + 2f 2 f 5 T u xT v + 2/ 2 / 6 7 n xT VI + 2/ 5 / 6 T v xT VI (Al) 



ff* — ► zz 

I. s-channel if annihilation 

II. s-channel /i annihilation 

III. t-channel f exchange 

IV. u-channel f exchange 

V. point interaction 



—92'mz cos(/3 — a)/ cosfl^Xfl^zsin 2 ^ cos(o; + cosd w ) 
—g-imz sin(/3 — a)/ cos w )(—g2m z sin 2 ^ sin(o; + /?)/ cos W ) 
-5(2 sin 2 ^/ cos# w ) 2 
-g 2 sin 2 6 w / cos9 w ) 2 
-2g\ sm A Q w / cos 2 9 w 

I - 4m|s + s 2 )/(4m|(m^ - s) 2 ) 
4m| S + S 2 )/(4m|(m 2 - S ) 2 ) 
m?„ - 4ml m| + 6m 4 . m| - 4m? m| + m| - 4m? t + Amfm 2 z t + 



Tfl Tfl Z 1 T/J Z Ttf Z ' Z TR T R 

Am 2 fR m z t - Am%t + 6m^t 2 + Am 2 fR m 2 z t 2 + 6m|t 2 - Am\ R t 3 - Am 2 z t 3 + t 4 ) / 
m z (m ffl 



12ml 



12m| 



m f R 



Aml^m z + 6m~ m| - 4m 2 m| + m| - Am%„u + Amf n m 2 7 u + 



Am 2 fR m z u - Am%u + 6m~ fl w 2 + Am 2 R m 2 z u 2 + 6m|w 2 - Am 2 R u 3 - Am 2 z u 3 + w 4 )/ 

12m| -4m|s + s 2 )/(4m|) 

12m| - 4m|s + s 2 )/(4m|(m^ - s)(mj; - s)) 

—&m TR m 2 z — 20m 2 fl m| — 6m| + rn\ R s + 2m\ R m 2 z s + 5m|s + 8m 2 R m 2 z t + 
8m|t - 2m 2 R st - 2m 2 z st - 2m 2 z t 2 + st 2 + Am\ R m 2 z u + 4m|w - Am 2 z tu) / 
2m z (m 2 H -s)(m 2 R -t)) 

—§m\ R m 2 z — 2Qm 2 R m z — 6m| + rnf R s + 2m 2 fl m|s + 5m|s + Am 2 R m 2 z t + 

/I 00 /I O O O 00 0\ / 

Am z t + 8m fR m z u + 8m z u — 2m~ R su — 2m z su — Am z tu — 2m z u + su )/ 
2m z (m 2 H - s)(m 2 _ - «)) 



12m| - 4m|s + s 2 )/(Am z (m 2 H - s)) 
-6m~ R m 2 z - 20m 2 fR m z 



Qm 7 + m-s + 2m-m 7 s + 5m 7 s + 8m-m 7 t + 
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8m%t - 2m\ R st - 2m 2 z st - 2m 2 z t 2 + st 2 + Am\ R m 2 z u + Am z u - Am 2 z tu) / 
(2m z (ml-s)(m 2 R -t)) 
T u xTiy = (— 6mj R m 2 z — 2tim 2 R m z — 6m| + mj R s + 2m 2 R m 2 z s + 5m z s + Am 2 R m 2 z t + 

A A O O 00 0\ / 

Am z t + 8m fR m z u + 8m z u — 2m fR su — 2m z su — Am z tu — 2m z u + su )/ 
{2m\{m\ - s)(m 2 R -u)) 
T u xT Y = (12m z - Am 2 z s + s 2 )/(Am z (m 2 h - s)) 
T ul xT w = {m% R + Y2m% R m 2 z + ?>8m\ R m\ + \2m\ R m\ + m| - Am\ R m 2 z s - 2Am 2 R m\s - 
Am%s + Am z s 2 - 2m% R t - Umj R m 2 z t - Um 2 fR m z t - 2m\t + Am\ R m 2 z st + 
Am z st + mj R t 2 + 2m 2 R m 2 z t 2 + m z t 2 - 2m% R u - \Am\ R m 2 z u - lAm 2 R m%u - 

f\ 00 A A 00 A 

2m z u + Am fR m z su + Am z su + Am- R tu + lQm fR m z tu + Am z tu — Am z stu — 
2m\J 2 u — 2m 2 z t 2 u + m\ R u 2 + 2m\ R m 2 z u 2 + m^u 2 — 2m 2 R tu 2 — 2m 2 z tu 2 + 
t 2 u 2 )/(m z (m 2 R -t)(m 2 R -u)) 

- / /10 0/1 00 /I 00 

T m xTv = (— 6m fR m z — 20m fR m z — 6m z + m fR s + 2m~ R m z s + 5m z s + 8m fR m z t + 
8m|t - 2m 2 R st - 2m|st - 2m|t 2 + st 2 + Am\ R m 2 z u + 4m|w - Am 2 z tu) / 
(2m z (m 2 R -t)) 

TiyxTy = (— Qmf R m 2 z — 2Qm 2 R m z — 6m| + w4 fl s + 2m 2 R m 2 z s + 5m^s + Am 2 R m 2 z t + 
Am z t + 8m\ R m 2 z u + 8m z u — 2m 2 ~ R su — 2m 2 z su — Am 2 z tu — 2m 2 z u 2 + su 2 )/ 
(2m z {m 2 R -u)) 

\T\ 2 = f 2 T 1 xT 1 + f 2 T u xT u + flT m xT m + f 2 r w xT w + f 2 T v xT Y + 2f^^^^ 

2f l hT l xT lll + 2}\}\T l xT lY + 2/ 1 / 5 7JxT v + 2f 2 f 3 T u xT m + 2/ 2 / 4 7M v + 
2f 2 f 5 T u xT v + 2f 3 f 4 T m xT w + 2f 3 f 5 T m xT v + 2f 4 f 5 T w xT v (A2) 

f f * ► Z7 

I. t-channel f exchange 

II. u-channel f exchange 

III. point interaction 



h 


= e(—g 2 sin 2 9 w / cos9 w ) 






= e(— g 2 s'm 2 9 w / cos9 w ) 




h 


= 2eg 2 sin 2 9 w / cos 9 W 




Tjx% 


= (-2m 6 fR + Am\m 2 z - 2m 2 m z + 2m\ R t + 8m 2 m 2 z t - 
Am 2 z t 2 -2t")/(m 2 z (m 2 R -t) 2 ) 


- 2m z t + 2m 2 ~ R t 2 + 


TuxTu 


= (~ 2m f R + Am i R m z - 2m 2 ~ R m z + 2m\ R u - 8m 2 R m 2 z u - 
Am 2 z u 2 -2u 3 )/(m 2 z (m 2 R -u) 2 ) 


- 2m z u + 2m\ R u 2 + 


r m xT m 


= 3 
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TixTn = (&n4 R + 36mf R m 2 z + 6m 2 fi m| - 2m% R s - 24m 2 fR m 2 z s - 2m z s + Am 2 z s 2 - 
7mj R t - Ylm\ R m 2 z t - m z t + 2m\ R st + 4m 2 z st + m\ R t 2 + m 2 z t 2 - 7m\ R u - 
12mf — m z u + 2m 2 R su + Am 2 z su + Sm 2 R tu + 2m 2 z tu — 2stu — t 2 u 
+m 2 fR u 2 + m 2 z u 2 - tu 2 ) I \m 2 z (m 2 R - t)(m 2 - R - u)) 

TjxTui = {—2m\ R — lbm 2 R m 2 z — 3m z + m 2 fR s + 5m 2 z s + 2m 2 ~ R t + 3m 2 z t — st + 
2m\ R u + 4m|w - 2tu) / {2m 2 z (m 2 ~ R - t)) 

T n xT m = {-2m\ R - \hm\ R m\ - 3m| + m\ R s + 5m|s + 2m 2 R t + Am 2 z t + 2m\ R u + 
3m 2 z u — su — 2tu) I (2m 2 z (m 2 R — u)) 
\T\ 2 = / 1 2 ^x^ + / 2 2 7i I x^ I + /3 2 7i„x^„ + 2/ 1 / 2 ^x7i I + 2/ 1 / 3 7ix^ II + 

2f 2 f3T u xT m (A3) 

ff* — > 77 

I. t-channel f exchange 

II. u-channel f exchange 

III. point interaction 

h = e 2 
h = e 2 
h = ~2e 2 
%x% = (4mt R + 8m 2 R t + 4t 2 )/(m 2 R -t) 2 
TnxTn = (4mj R + 8m 2 R u + 4u 2 )/(m 2 R -u) 2 
T m xT m = 4 

TjxTn = (36mf R - 24m 2 fR s + 4s 2 - V2m\ R t + Ast + t 2 - \2m\ R u + Asu + 2tu + u 2 ) / 

((m 2 R -t)(m 2 R -u)) 
%xT m = (-12m 2 R + 5s + 4u)/(2(m 2 R -t)) 
T u xT m = (-12m 2 R + 5s + At)/(2(m 2 R -u)) 

\T\ 2 = f 2 T 1 xT I + flT u xT u + f 2 T m xT m + 2f 1 f 2 T 1 xT n + 2f 1 f 3 T 1 xT m 

2f 2 fsT u xT m (A4) 

ff* — > Zh[H] 

I. t-channel f exchange 

II. u-channel f exchange 

III. s-channel Z annihilation 

fi = (— #2 sm 2 9 w / cos9 w )(— g2m z sin 2 ^ sin[— cos] (a + (5)/ cos6 w ) 
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h = 


— {—92 sm 2 9 w / cos6 w )(—g2m z sin 2 ^ sin[— cos] (a + 13)/ cos9 w ) 


fs = 


(—g 2 sin 2 9 w / cos9 w )(—g 2 m z sin[cos](/3 — a)/ cos 9 W ) 


r Y x% = 


(mj R + (m| - t) 2 - 2m 2 R {m 2 z + t))/{m 2 z {m 2 R - t) 2 ) 


TuxTu = 


(m 4 fR + (m| - u) 2 - 2m 2 fi (m| + u))/{m 2 z {m 2 R - u) 2 ) 


TixTn = 


(mf R + m z + rn 2 fR {Qm 2 z — t — u) + tu — m 2 z {2s + t + u))/ 




(m 2 z (m 2 R -t)(m 2 R -u)) 


TjxTm = 


(t(t -u) + m 2 fR (-8m 2 z -t + u) + m 2 z (2s -t + u))/ 




(2m 2 z (m 2 z -s)(m 2 fR -t)) 


T u xT m = 


((t - u)u + m 2 fR (8m 2 z - t + u) + m 2 z (-2s -t + u))/ 




(2mi(mi -s)(m 2 fR - u)) 


T m xT m = 


(-lQm 2 fR m 2 z + Am 2 z s + (t - u) 2 )/(4m 2 z (m 2 z - s) 2 ) 


\T\ 2 = 


fl%x% + f 2 T u xT u + flT m xT m + 2f 1 f 2 T 1 xT u + 2f 1 f 3 T 1 xT UI + 



2f 2 f 3 T u xT m 

ff* — ► <yh[H] 

I. t-channel f exchange 

II. u-channel f exchange 



fi = (e)(— g 2 m z sin 2 ^ sin[— cos] (a + (3)/ cos9 w ) 

f 2 = —(e)(—g 2 m z sin 2 ^ sin[— cos] (a + (3)/cos9 w ) 

%x% = -2{m 2 R +t)/{m 2 R -t) 2 

TjxTu = -(-6m 2 R + 2s + t + u)/((m 2 R -t)(m 2 R -u)) 

T u xT u = -2(m 2 R +u)/(m 2 R -u) 2 

\T\ 2 = f 2 T l xT 1 + f 2 T u xT u + 2f 1 f 2 T 1 xT u 



ff* > ZA 



I. s-channel h exchange 

II. s-channel H exchange 



fi = (g 2 cos(a — /3)/(2cos9 w )(— g 2 m z sin 2 9 w sm(a + (3)/ cos9 w ) 

/2 = (g 2 sin(a — (3)/ (2 cos 9 w )(g 2 m z sin 2 9 w cos(a + (3)/ cos 9 W ) 

T Y x% = (m A + (m 2 z - S ) 2 -2m 2 A (m 2 z + S ))/(m 2 z (m 2 -s) 2 ) 

TnxTn = (m\ + (m| - s) 2 - 2m 2 A (m 2 z + s))/(m 2 z (m 2 H - s) 2 ) 

TjxTn = (m\+ (m 2 z -s) 2 -2m 2 A (m 2 z + s)/(m 2 z (m 2 H -s)(m 2 h -s))) 

\T\ 2 = f 2 T 1 xT l + f 2 r ll xT 11 + 2f 1 f 2 T l xT ll 
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TT > TT 

III. s-channel Z annihilation 

IV. s-channel 7 annihilation 

V. t-channel x exchange 



f?,c 


= ( -g 2 sin 2 9 w j cos 9 W ) (g 2 {l - 4sin 2 9 w ) /(4cos9 w )) 


hd 


= (-g 2 sm 2 9 w / cos 9 w )(-g 2 / (4 cos0 w )) 


he 


= e 2 


K 


= gi Nn/V2 


K' 


= giN n /V2 




- o(-4 f 2 m 2 s - 4 f 2 m 2 s f 2 s 2 4- f 2 s 2 - f 2 t 2 - f 2 / 2 + 2 f 2 tv + 2 f 2 hi 




flu 2 - f 2 3d u 2 )/(m 2 z - s) 2 


T m xT lv 


= 2hcfiMm 2 R s -s 2 + t 2 - 2tu + u 2 )/((m 2 z - s)s) 




= 2f 2 c (-4m 2 R s + s 2 -t 2 + 2tu - u 2 )/s 2 


r m xT v 


= ~2(hc + hd)K 2 (4m 2 R s -s 2 + (t- «) 2 )/((m| - s)(m£ - t)) 




= 2(h c )K 2 (4m 2 R s -s 2 + (t- u) 2 )/(s(m£ - t)) 


TyXTy 


= (16K 2 K' 2 (mt R - tu))/{{m£ - t){-m^ 2 + t)) 


\T\ 2 


= T m xT m + T w xT iy + 2T m xT iy + 




4 

£ ((TixTy + TnxTy + T m xT y + T iy xT y )/2 + T v xT v ) 



f f * — > // 

III. s-channel Z annihilation 

IV. s-channel 7 annihilation 



/3c = (-g2sm 2 9 w /cos9 w )(g 2 (-2Ti + 4Q f sm 2 9 w )/(4cos9 w )) 

hd = (-g 2 sm 2 9 w /cos9 w )(g 2 (2Ti)/(4cos9 w )) 

he = -e/e 2 

T m xT m = 2(16f! d m 2 R m 2 f -4flm 2 R s-4f 2 d m 2 R ^ 

flf - fit 2 + 2/ 3 > + 2/ 3 > - f 2 c u 2 - flu 2 )/(m 2 z - s) 2 

T m xT w = 2hef^m 2 R s-s 2 + t 2 -2tu + u 2 )/{{m 2 z - s)s) 

T iy xT iy = 2f 2 c (-4m 2 fR s + s 2 - t 2 + 2tu - u 2 ) / s 2 

\T\ 2 = (T m xT m + T iy xT iy + 2T m xT iy )(x3 for quarks) (A9) 
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ff* — ► tt 

I . s-channel H annihilation 

II. s-channel h annihilation 

III. s-channel Z annihilation 

IV. s-channel 7 annihilation 



fia = (g 2 m z sin 2 9 w cos(a + (3)/ cos 9 w )(-g 2 m t sin a/ (2m w sin /3)) 

f2a = (-g2'mzsin 2 9 w sin(a + (3)/ cos9 w )(-g 2 rntcosa/(2m w sin/3)) 

he = (-g2sin 2 9 w / cos9 w )(g 2 (-l + 4Q t sin 2 9 w ) / (4cos9 w )) 

fu = (- g 2 sin 2 9 W / cos 9 w )(g 2 / (4 cos9 w )) 

fic = ~e t e 2 

7 Y x% = 2f 2 a (-4m 2 + s)/(-m 2 H + s) 2 

r n xT n = 2f 2 2a {-4m 2 t + s)/(-m 2 + s) 2 

T m xT m = 2(16/ 3 >?X - ±fWf R s ~ ±fW fR s - 4f 2 d m 2 s + f 2 c s 2 + f* d s 2 - 
fsf - flf + Vltu + 2/ 3 > - f 2 3c u 2 - f 2 3d u 2 )/(m 2 z - s) 2 

TtxTn = 2f la f 2a (-4m 2 + s)/((-m 2 H + s)(-m 2 h + s)) 

TixT m = 4f la f 3c m t (t - u)/((m 2 H - s)(-m 2 z + s)) 

T u xT m = 4f 2a f 3c m t (t-u)/((m 2 h - s)(-m 2 z + s)) 

TjxTiy = 4f la f Ac m t (t - u)/((m 2 H - s)s) 

TnxTiy = 4f 2a f Ac m t (t-u)/{{m 2 h - s)s) 

T m xT w = 2f 3c f Ac {4m 2 fR s - s 2 + t 2 - 2tu + u 2 ) / {{m 2 z - s)s) 

T w xT lv = 2f 2 c {-4m 2 ~ R s + s 2 -t 2 + 2tu-u 2 )/s 2 

\T\ 2 = 3 (TjxTi + T n xT u + T m xT m + T w xT w + 2%x% Y + 2%xT m + 2%xT lw + 

2T u xT m + 2T u xT lv + 2T m xT lv ) (A10) 



ff* — ► hh 

I. s-channel h annihilation 

II. s-channel H annihilation 

III. point interaction 

IV. t-channel f exchange 

V. u-channel f exchange 



/1 = (—g 2 m z sin 2 9 w sin(a + ft) / cos9 w )(—3g 2 m z cos2asin(a + (3) / (2cos9 w )) 
f 2 = (g 2 m z sin 2 9 w cos(a + (3)/ cos 9 w )(g 2 m z (cos2a cos(a + f3) — 

2 sin(2a) sin(a + (3))/ (2 cos 9 W )) 
f 3 = —g\ cos 2a sin 2 9 w /(2 cos 2 9 w ) 
fi = (— g2^z sin 2 9 w sin(o; + (3)/ cos9 w ) 2 
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h 




(— q 2 m z sin 2 # w sinfo; + 3)1 cos 9 W ) 2 


TixTi 


— 


(m\ - s)- 2 


TnxTn 




(m 2 H - s)- 2 


T m xT m 




1 


T w xT w 


= 


(m 2 R -t)- 2 


T v xT Y 


= 


(m 2 R -u)- 2 


TixT n 


— 


l/((m 2 H -s)(m 2 h -s)) 


TixT m 


= 


l/(m 2 h - s) 


TixT w 


= 


l/((m 2 h -s)(m 2 fR -t)) 


%xT v 


= 


1 /(( m l- s)(m 2 R -u)) 


r u xT m 


= 


l/(m 2 H - s) 


T n xT w 


= 


l/((m 2 H -s)(m 2 fR -t)) 


r u xT Y 






T m xT w 


= 


1/(^4 - 1) 


T ul xT v 




V( m L - u ) 

1 v 1 R 7 


T w xT v 




l/((m 2 fR -t)(m 2 fR -u)) 


\T\ 2 




f^xT, + f 2 T n xT n + f 2 T m xT m + f 2 T w xT w + f 2 5 T y xT v + 2f l f 2 T l xT Il + 
2f 1 f 3 T I xT m + 2f 1 }\T l x%j + 2f 1 f 5 TixT Y + 2f 2 f 3 T n xT m + 2f 2 f 4 T u xT w + 
2f 2 f 5 T u xT v + 2f 3 f 4 T m xT w + 2f 3 f 5 T m xT v + 2f 4 f 5 T w xT v (A) 



ft* — ► hA[HA] 

I. s-channel Z annihilation 



fi = (-92sin 2 6 w / cos6 w )(g 2 cos[sin](a - (3)/(2cos9 w )) 
T.xT, = (t-u) 2 /(m 2 z -s) 2 
\T\ 2 = flTixTj (A12) 



f f * ► W+H 



I. s-channel H annihilation 

II. s-channel h annihilation 



fi = (g 2 m z sin 2 6 w cos(a + (3)/ cos6 w )(—g 2 sin(a - (3)/2) 

/2 = (-g 2 m z sin 2 6 w sin(a + (3)/ cos6 w )(-g 2 cos(a - /3)/2) 

Tix% = K+ + K - s ) 2 - H+K + s ))/WK - s) 2 ) 

TjxTn = (ro 4 H+ + (m^ - s) 2 - 2mj, + (m 2 ff + s))/ (m^mjj - s)(m^ - s))) 
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T n xT n = {m 4 H+ + (m 2 w - sf - 2m 2 H+ (m 2 w + s))/{m 2 w {m 2 h - s) 2 ) 

\T\ 2 = f 2 T l xT 1 + f 2 T l xT I + 2f 1 f 2 T 1 xT n (A13) 



f f * ► AA 

I. s-channel if annihilation 

II. s-channel h annihilation 

III. point interaction 



fi = (g 2 mz sin 2 9 w cos(a + (3)/ cos9 w )(g 2 mz cos2(3cos((3 + a)/(2cos6 w ) 

f2 = (—g 2 m z sm 2 9 w sin(o; + P)/ cos 9 w )(—g 2 m z cos2/3sin(/3 + a)/(2cos9 w ) 

f 3 = —g 2 cos 2/3 sin 2 9 W / (2 cos 2 9 W ) 

TixTi = (m 2 H -s)- 2 

TnxTn = (m 2 - s)' 2 

T m xT m = 1 

T.xTn = \/{{m 2 H -s){m 2 h -s)) 

TjxTm = \/{m 2 H -s) 

T n xT m = l/(m 2 h -s) 

\T\ 2 = f 2 T l xT 1 + f 2 TuxT u + fiT ul xT m + 2f 1 f 2 T l xT u + 2f 1 f 3 T l xT m + 

2f 2 hT ll xT m (A14) 



ff* — > hH 

I. s-channel H annihilation 

II. s-channel h annihilation 

III. point interaction 

IV. t-channel r exchange 

V. u-channel f exchange 



fi = (g 2 mz sin 2 9 w cos(a + ft)/ cos 9 w )(g 2 mz(2sm2a cos((3 + a) + 

sin(/3 + a) cos2a)/(2 cos9 w )) 
f2 = (—g 2 m z sm 2 9 w sm(a + (3)/cos9 w )(—g 2 mz(2sm2asm((3 + a) — 

cos(/3 + a) cos 2a) / (2 cos 9 W )) 
fs = — g 2 sin 2a sin 2 9 W /(A cos 2 9 W ) 

/ 4 = (^2"^z sin 2 ^ cos(o; + (3)/ cos9 w )(— g 2 m z sin 2 ^ sin(o; + (3)/ cos9 v 
h = {g 2 m z s\n 2 9 w cos(a + (3)/ cos9 w )(-g 2 m z sm 2 9 w sin(a + (3)/ cos9 v 
%x% = ( m 2 H -s)- 2 

T.xTu = l/((m 2 H -s)(m 2 -s)) 

TjxTm = \/{m 2 H -s) 
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%xT lY 


= l/((m 2 H -s)(m 2 fR -t)) 




= l/((m 2 -s)(m 2 fR -u)) 


TuxTu 


= (m 2 h -s)- 2 

\ it / 


TuxTm 


= l/(™ 2 h -s) 


T n xT w 


= l/((m 2 h -s)(m 2 R -t)) 




= l/((m 2 H -s)(m 2 fR -u)) 


T m xT m 


= 1 


T m xT lY 


= v«-o 


T m xTy 


= l/(m 2 R -u) 


T 1Y xT lY 


= (ml R -t)- 2 


T lY xT Y 


= l/((m 2 fR -t)(m 2 fR -u)) 


TyXTy 


= (ml R -u)- 2 


\T\ 2 


= fftxTj + flTnxTu + f 2 T m xT m + f 2 A T lY xT lY + f 2 T Y xT Y + 2hf 2 %xT n + 



2f 1 f 3 T l xT m + 2f 1 f 4 T 1 xT lY + 2f 1 f 5 T l xT Y + 2f 2 f 3 T u xT m + 2f 2 UT n xT 1Y + 
2f 2 f 5 T u xT Y + 2f 3 f 4 T m xT lY + 2f 3 f 5 T m xT Y + 2f 4 f 5 T w xT Y (A15) 

ff* — ► HH 

I. s-channel H annihilation 

II. s-channel h annihilation 

III. point interaction 

IV. t-channel f exchange 

V. u-channel f exchange 



fi = (d2 m z sin 2 ^ cos(a + (3)/ cos9 w )(— 3g 2 m z cos 2a cos(a + (3)/ (2 cos9 w )) 

h = (— 92'mz sin 2 ^ sin(o; + (3)/ cos 6 w )(g 2 m z (cos 2a sin(o; + (3) + 
2 sin 2a cos(a; + (3))/ (2 cos 6 W )) 

fs — glcos2asm 2 6 w /(2cos 2 9 w ) 

/ 4 = (g 2 m z sm 2 6 w cos(a + (3)/cos9 w ) 2 

/s = (g2'TTiz sin 2 6 W cos(a + (3)/ cos 6 W ) 2 
|T| 2 = / 1 2 7ix7i + / 2 7I I x7I I + /3 2 ^ n x7I n + / 4 2 7ivx7iv + / 5 2 T v xT v + 2/ 1 / 2 7Ix7i I + 
2f 1 f 3 T l xT lll + 2f l f 4 T l xT lY + 2f 1 f 5 T l xT Y + 2f 2 f 3 T u xT m + 2f 2 f 4 T u xT IY + 
2f 2 f 5 T n xT Y + 2f 3 UT m xT lY + 2f 3 f 5 T m xT Y + 2f 4 f 5 T w xT Y (A16) 

The 7[x7[ . . . are the same as for ff* — ► hh, with (m^ <-> m#). 
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ff* — ► H+H 



I. s-channel annihilation 

II. s-channel h annihilation 

III. s-channel Z annihilation 

IV. s-channel 7 annihilation 

V. point interaction 



h 


— 


(q 2 m z sm 2 9 w cosfa + 0)1 cos 6 w )(—q 2 (mw cos(/3 — a) — 
m? cos 23 cos(3 + a) 1(2 cos 6 W )) ) 


h 


— 


(—92'mz sin 2 9 w sinfo; + (3)1 cos 9 W ) (—g 2 (mw sin(3 — a) + 
mz cos2/3sin(/5 + a)/(2cos 9 W ))) 


h 


— 


(-g 2 sm 2 9 w / cos 9 w )(-g 2 cos 2^/ (2 cos 


U 


— 


-e 2 


h 


— 


—ql cos 23 sm 2 9 w l (2 cos 2 6> w ) 




— 


(m 2 H - s)- 2 


TixTn 


— 


l/((m 2 H -s)(m 2 h -s)) 

1 \ \ ±1 / \ lb / / 


TixTm 


— 


(t - u)/((m 2 H - s)(m| - s)) 




— 


(t-u)/(-(m 2 H s) + s 2 ) 


TlXTy 


= 


l/(m 2 H - s) 


T n xT n 




(m 2 h - s)- 2 


TiixTui 




(t - u)/((m 2 h - s)(m 2 7 - s)) 


T u xT 1Y 




(t-u)/(-(m 2 h s) + s 2 ) 


T u xTy 




l/(m 2 - s) 


T m xT m 




{t-uf/{m 2 z - s) 2 


T m xT w 




-{{t-u) 2 /{{m 2 z -s)s)) 


T m xTy 




(t-u)/(m 2 z - s) 


T w xT w 




(t - uf/s 2 


T w xT y 




(-t + u)/s 


%jxT w 




1 


\T\ 2 




fl%x% + f 2 T u xT u + fiT m xT m + flT lw xT lw + / 2 T v xT v + 2/ 1 / 2 7 I x7 II + 
2hh%xT m + 2/ 1 / 4 7ix7rv + Vih%xT y + 2f 2 f 3 T u xT m + 2f 2 f 4 T n xT w + 
2/ 2 /5^ixT v + 2f 3 UT m xT w + 2f 3 f 5 T m xT v + 2f 4 f 5 T w xT v (AJ 



rr — ► rr 

I. t-channel x exchange 

II. u-channel x exchange 
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K = gi N n /V2 
K' = g 1 N jl /y/2 

l xi m xj s )/(( m x 2 i - u )( m x 2 j 
7JxT„ = (16K*K' 2 m^m^s)/((mi-t)(mt 2 

4 



TixTj = (16K 2 K' 2 m^m^s)/((m^ -t)( mi 2 
T n xT„ = (lQK 2 K' 2 m^m^s)/((m^ -u)(m^ 2 



\T\ 2 = £(Tp<T^T„xT„ + 27JxT„) 



ft — > t£ 

I. t-channel x exchange 



K = 9l N n /V2 

K' = g 1 N jl /V2 

TixTi = -16K 2 K' 2 (m 2 R m 2 iR -tu)/((m^ 2 -t)(m 

\T\ 2 = £71x71 



fl — > t£ 

I. t-channel x exchange 



K = g 1 N il /V2 
K' = g 1 N jl /V2 

4 



|7f = J2 T ^ T i 



fx — ► 



I. s-channel r annihilation 

II. t-channel f exchange 



/! = -(^iN jl /V2)(- 52 /(2cosM) 
/ 2 = (^i/v^)(-^siii 2 w /cos0 w ) 
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TixTi = 2(2 sm 2 O w ) 2 {m\s - mf R s + s(—m% + m%s + m z t - st - m 2 z u) + 

m 2 R (2m z — 2m 2 z s + s 2 + st + su) — m^(2m z — 2m 2 z s + s(t + u))) / [m 2 z s 2 ) 
TnxTn = 2{m\ - t){m\ R + (m| - t) 2 - 2m 2 fR (m 2 z + t))/(m 2 z (m 2 fR - t) 2 ) 
TixTn = -{2sin 2 6 w ){ml R {s + t-u)+m\{-{m 2 z s)+m 2 z t + st - t 2 + 

m 2 R (8m 2 z — s + t — u) — 5m 2 z u + tu) + (m| — t){m 2 z s — s 2 + m z t — t 2 — 
m 2 z u + u 2 ) + m\ R {2m 2 z s - s 2 - 2m 2 z t - st - 2t 2 - 2m 2 z u + tu + u 2 ))/ 
(m 2 z s(m 2 fR - t)) 

\T\ 2 = f 2 T 1 xT 1 + f 2 T u xT ll + 2f 1 f 2 T 1 xT ll (A21) 
TX ► 

I. s-channel r annihilation 

II. t-channel f exchange 

h = -{9iN 3l /V2){e) 

T^Ti = A{m\-m\ R - su-m\{t + u)+m 2 fR {s + t + u))/s 2 

r n xT n = -^m\-t){m\ R + t)/{m 2 fR -tf 

TixTn = (s 2 + t 2 -u 2 + ml R (-s + 3t + u) + ml(-8ml R + s-t + 5u))/( y s( y ml R -t)) 

\T\ 2 = f 2 T 1 xT l + f 2 T u xT ll + 2f 1 f 2 T l xT ll (A22) 

fx ► r/i[if] 

II. t-channel f exchange 

/2 = 

T u xT u = 
\T\ 2 = 
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